Search for ?7-mesic nuclei in recoil-free transfer reaction 
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Abstract 

We have studied the reaction p + '^'^Al ^He + p + tt~ + X at recoil-free kinematics. An 77 
meson possibly produced in this reaction would be thus almost at rest in the laboratory system and 
could therefore be bound with high probability, if nuclear r] states exist. The decay of such a state 
through the A^*(1535) resonance would lead to a proton- vr" pair emitted in opposite directions. 
For these conditions we find some indication of such a bound state. An upper limit of w 0.5 nb is 
found. 
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The study of A and S hypernuclei, which are nuclear bound systems of short-hved 
hadrons, has proved to be a very useful tool for gaining information about A — and 
T, — N interactions in nuclei. Also, studies of 7r~- and K~- atomic levels have provided 
useful information about tt — N and K — N interactions in nuclei. However, there has never 
been, to the best of our knowledge, an observation of a neutral pseudoscalar meson bound 
strongly in the nucleus. Observation of such bound states would open new possibilities in 
nuclear and particle physics with respect to the structure of such nuclei, the rjNN* coupling 
constant and the behavior of the Sn nucleon resonance in nuclei. 

In contrast to the pion-nucleon interaction, the //-nucleon interaction at small momenta 
is attractive and sufficiently strong. This attraction can be seen from the fact that the rj 
threshold (1488 MeV) is situated just below the A^*(1535) resonance which couples strongly 
to the rj — N channel. Initial calculations by Bhalerao and Liu [l| obtained attractive s-wave 
r] — N scattering lengths o^at = (0.28 + 0.19i) fm and a^^N = (0.27 + 0.22i) fm, using the 
Tx — N phase shifts calculated by Arndt and the CERN theory group, respectively. With 
these phase shifts, Haider and Liu l2l have shown that r? can be bound in nuclei with A > 10. 

n fi ri 

Other groups have also found similar results [3|, U, 15|. Recent analyses of the experimental 
data and different theoretical calculations predict a range of values for the rj — N scattering 
length from 0.2 fm to 1.0 fm for the real part, and from 0.2 fm to 0.35 fm for the imaginary 
part. The higher values for the real part of o^at have led to speculations that the //-bound 
state might be possible even for lighter nuclei. An overview of this topic is given in Ref. Q]. 
There have been previous searches for the proposed ?7-mesic nucleus. First experiments 

n n 

searching for //-mesic nuclei at BNL [7[ and LAMPF |8[ by using a missing-mass technique 
in the (vr"*",/)) reaction came to negative or inconclusive results. Later it became clear that 
the peaks are not necessarily narrow and that a better strategy of searching for ?7-nuclei is 
required. Furthermore, the BNL experiment was in a region far from the recoilless kinemat- 
ics, in which the cross section is substantially reduced 9|]. More recently, the existence of 
?7-mesic ^He was claimed to have been observed in the reaction y^He n^pX using the 

H)ton beam at MAMI 1^. It has, however, been pointed out 11| that the data of Ref. 
does not permit an unambiguous determination of the existence of a '^ife?7-bound state. 
The suggestion that ^Hei] is not bound is also supported by the theoretical studies of Refs. 
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The present experiment makes use of the transfer reaction 

p + ^Z ^"^He + ^-^Z -1)^7] (1) 

The ^He were measured at zero degrees with the magnetic spectrograph Big Karl [l^ by ray 
tracing in the focal plane with two packs of multi-wire chambers. This was followed by two 
hodoscope layers, separated by 4 m which provided an additional time-of- flight measurement. 
The beam momentum (pboam = 1745 MeV/c) and the setting of the spectrograph were chosen 
such that for binding energies in the range 0-20 MeV, the t] is produced almost at rest. The 
^He spectrum is expected to be dominated by particles being emitted during the nuclear 
cascade process. In order to reduce this background a coincidence was required among ^He 
and events produced through a second step 

rj + n —>■ n~ + p. (2) 

Because the overall rjn system (or A^*) is almost at rest, energy and momentum conservation 
require the two charged particles to be emitted back-to-back to each other with energies of 
!^100 MeV for the proton and ^^^348 MeV for the pion. Such a clear pattern is smeared out 
by Fermi motion resulting in a distribution around ^150 ° with a width of 40°. For the 
measurement of these particles a dedicated detector ENSTAR was built, details of which are 
described in Ref. jl^. Briefly, it consists of three cylindrical layers of scintillating material 
surrounding the target. Each layer is divided into long bars thus allowing a measurement 
of the azimuthal emission angle. The bars of the middle layer are further divided along 
the length in order to measure the polar emission angle. While the protons of interest are 
stopped in the middle layer of the detector, pions pass through all layers giving only AE 
information. 

Although, some calculations predict ^He to be large enough to bind rj mesons, Garcia- 
Recio et al. expect more medium mass nuclei (A ~ 24) to show stronger binding. On 
the other hand, heavier nuclei will have broader states making them harder to detect on 
a smooth background. Furthermore, the final nucleus should not have too many excited 
states, which is the case for even-even nuclei. The ideal target should thus be odd-odd, but 
such a nucleus does not exist as a solid target so we were limited to an odd-even system. 
As a compromise among these different factors we choose ^'^Al. The target thickness of 
1 mm, corresponding to a resolution of 2 MeV, was chosen in order not to spoil the natural 
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width of the bound state. Two runs were performed with different spectrometer momentum 
settings (po = 859 MeV/c and 897 MeV/c). An integrated luminosity of 0.50 ± 0.05 pb"^ 
was accumulated for each run. 

Prior to the experiment, the ENSTAR detector was calibrated as described in Ref. jl^ . 
Due to the high brilliance proton beam the experiment was performed with minimal back- 
ground even though the ^He were measured in the forward direction. 
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FIG. 1: Energy loss AE as function of the time-of-flight (TOF) for particles in the focal plane area 
of the magnetic spectrograph. 



In Fig. [T] the energy loss in the first hodoscope layer is shown as a function of the time- 
of-fiight. Different particle groups can be clearly identified. Beam particles do not enter the 
focal plane because their charges differ by a factor of two. This would not be the case in 
a deuteron induced reaction where break-up protons would flood the focal plane detectors. 
The inclusive ^He spectra are uniformly distributed when the data are corrected for the 
acceptance of the spectrograph. ^He within the angles ~ 3° in the vertical and ~ 0.6° in 
the horizontal direction were recorded by the spectrograph focal plane detectors. 

The coincidence required between the focal plane detectors and the ENSTAR detector was 
achieved by measuring the time between the flrst hodoscope layer and one of the individual 
ENSTAR elements. A peak-to-background ratio of 3.2:1 was obtained and the background 
was subtracted. 

For the beam momenta used, the selection of ^He means that the residual system is at rest 
with an excitation energy of ~ 550 MeV. The only background which could give the same 
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pattern as the N* decay would be a deuteron, stopping in the middle layer in association 
with a higher energy proton punching through all detectors. A gate was therefore put on 
pions on a AE — E spectrum for events going through all layers. With such geometrical 
selections we obtain the missing-mass spectra for the two spectrograph settings shown in 
Fig. O The counts have been corrected for the acceptances of the spectrograph for the two 
settings. In order to minimize systematical uncertainties in areas of small acceptance, only 
the regions with acceptance 4% around the central momentum value have been retained. 
This eliminated data in less than 5% of the missing mass range. Applying the cuts from 
ENSTAR corresponding to 77 + n — * tt" + p leads to reduction in yield by a factor of ~ 10'^. 
Positive values of the binding energy BE correspond to the free or unbound rj production. 
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FIG. 2: (Color online) Missing mass spectra for two spectrograph settings indicated in the fig- 
ure. The counts were generated from acceptance-corrected ^He spectra measured in the magnetic 
spectrograph with two charged particles detected in the ENSTAR detector, which show the decay 
pattern of an N*. The sohd Une indicates zero excitation energy of a '^^Mg^ij system, i. e. binding 
is to the left of the Une. 
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FIG. 3: (Color online) Data in the peak region as function of the missing mass (upper abscissa) or 
the corresponding binding energy (lower abscissa). The upper panel shows the data without the 
back-to-back correlation. The solid line in the upper panel is a fitted constant to the whole data 
set. The solid curve in the lower panel is a constant fitted as background and gaussian on top of 
this background. The dotted curve is a fitted polynomial as background and a gaussian on top 
of this background. The insert shows the total data set for the back-to-back condition and two 
different background fits to the data outside the peak region. The errors are asymmetric due to 
the underlying Poisson statistics. 

Due to the large width of the A^* resonance of 100 to 200 MeV [l^, the yield should rise 
with phase space, while at the eta-mesic formation threshold it is expected to be zero. 
An indication of such a rise is somewhat seen in the data as demonstrated by the fitted 
polynomial shown in the insert of Fig. El which gives a better value than the one that is 
obtained by a fit with a constant. 

For both settings, there appears to be an enhancement from the threshold for "^^Mg ® ?7 
which is -23.8145 GeV to ^ —23.79 GeV. One may attribute all the counts to background. 
However, background should be randomly distributed and hence it is very unlikely for two 
different measurements to show the same structure. We, therefore, conclude, the structures 
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could be from the bound rj. 

In order to elucidate this point further we discuss the spectra in more detail. In Fig. [3] 
we show the binding energy spectra combined for the two settings. Since for both settings 
the same luminosity was acquired, the weighted arithmetic mean was used in overlap region. 
The figure shows the data without the strong back-to-back correlation requirement (upper 
panel) and with the requirement (lower panel). The unconstrained data do not show any 
structure and can be well described by a constant. For the data in the lower spectrum, 
the N* decay pattern is required. The counts are typically lower by an order of magnitude 
than that in the unconstrained case. Although, the N*(1535) can also decay with two-pion 
emission, this branch is small compared to the p7T~ channel 15|. 

The data show an enhancement around BE ^ —13 MeV. The significance of this structure 
is extracted according to the two methods given in Refs. [l6| and [? ] respectively. At first, 
we test the hypothesis of peak structure being fluctuation of background, i.e. the origin 
of the background is taken to be independent of the signal. The background outside the 
peak region, for simplicity approximated by a constant, was found to be 5.8±0.64. The 



significance [16|] is then given by (A^ — BG) / \JBG + asc where A^ is the total counts in the 
region of interest, BG is the total background in this region as determined from the fit to 
the outside region and aBc is error in the estimation of background value as taken from 
the fit. This yields a value of significance which is 5.3cr. Here we have assumed Gaussian 
errors. For the assumption of Poisson errors with asymmetric error bars (see Fig. [3]) the 
background is 6.2±1.0. This larger value is typical for Poisson distribution and hence the 
significance reduces to 4.9(7. Finally a Gaussian on top of the background was fitted to the 
whole data set. This yielded for the case of Poisson statistics 6.4±0.96 for the background, 
8.3±3.6 for the amplitude, -12.0±2.2 MeV for the centroid and 4.7±1.7 MeV for the width. 
The corresponding curves are shown in Fig. [HI Also a third order polynomial was fitted to 
the data. The resulting curves are also shown in Fig. [3l 

In the second method, the statistical significance is extracted by assuming the background 
events as well as the peak events on top of the background being Poisson distributed. Again a 
constant background and a Gaussian was assumed. A fit was performed using the maximum 



likelihood method. The significance is then defined as, ^/ —2AlnL. Here, AlnL is the 
difference in the values of the logarithm likelihood function with signal fixed to zero and at 
the best fit value. In this way, we obtain a value of 6.20a for the significance, assuming a 
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simultaneous determination of amplitude, centroid and width of the signal. The fit gives for 
the linear background 6.38 ± 0.53 together with the values, for the signal amplitude 8.55 ± 
3.05, for the centroid —13.13 ± 1.64 MeV and for the width 4.35 ± 1.27 MeV corresponding 
to a FWHM of 10.22±2.98 MeV. These results compare favorably with those from the first 
method. We, therefore, consider the present experimental results to provide a strong hint 
of a nuclear t] bound state. 

This allows us to give an upper bound for the cross section. With an estimated efficiency 
due to detector geometry and analysis selections of 0.70 ± 0.07 we find a = 0.152 ± 0.054 
(stat) ± 0.021 (syst) nb. If this "structure" corresponds to a bound rj decaying via reaction 
(Eqn. [2]), the cross section would be 0.46 ± 0.16 (stat) ± 0.06 (syst) nb, assuming an isospin 
branching ratio of 1/3. This cross section value can be compared with the "elementary" 



pd — > ^Het] reaction which, for the present 



3eam energy has a cross section integrated over 
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the spectrograph acceptance of ~ 39/i b 

In summary, we have measured the reaction p+'^'^ Al ^He + ir^ +p + X. The ^He ions, 
which were detected at zero degrees with a magnetic spectrograph, carried the beam momen- 
tum (recoilless kinematics). The remaining system has the mass m(^^Al) +m(?7) +BE with 
BE the binding energy. The ti~ + p system, measured with the ENSTAR detector, decays 
almost back-to-back with energies corresponding to an A^*(1535) at rest. The most probable 
scenario for the rj decay is through forming an in-medium A^* resonance which decays into p 
and 7r~. In this case, the remaining system is m(^^Al + m(A^*(1535)in medium)- Although in- 
stead of an ?7 a pion could be produced in the intermediate step forming another A^* nearby, 
which could lead to the back-to-back -K^p events. However, simple kinematical calculations 
show that this would require the momentum of a target nucleon to at least 210 MeV/c for 
which the probability is very low. In two spectra, taken at different spectrograph settings, an 
enhancement was found for negative binding energies close to the free production threshold. 
This is exactly what is expected from the bound rjN system. The enhancement may not be 
purely due to binding in the ground state only but also to an excited '^^Mg state. However, 
this requires pick up of more deeply lying nucleons which may be less likely than pick up 
of the least-bound nucleons. Binding energy spectra without strong A^* constraint do not 
show the enhancement. 
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